Mergers of black-hole binaries with aligned spins: Waveform characteristics 
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We conduct a descriptive analysis of the multipolar structure of gravitational-radiation waveforms 
from equal-mass aligned-spin mergers, following an approach first presented in the complementary 
context of nonspinning black holes of varying mass ratio [Baker et al. Phys. Rev. D 78, 044046 
(2008)]. We find that, as with the nonspinning mergers, the dominant waveform mode phases evolve 
together in lock-step through inspiral and merger, supporting the previous waveform description in 
terms of an adiabatically rigid rotator driving gravitational-wave emission - an implicit rotating 
source (IRS). We further apply the late-time merger-ringdown model for the rotational frequency 
introduced in Baker et al. (2008), along with an improved amplitude model appropriate for the 
dominant (2, ±2) modes. This provides a quantitative description of the merger-ringdown wave- 
forms, and suggests that the major features of these waveforms can be described with reference 
only to the intrinsic parameters associated with the state of the final black hole formed in the 
merger. We provide an explicit model for the merger-ringdown radiation, and demonstrate that this 
model agrees to fitting factors better than 95% with the original numerical waveforms for system 
masses above ~ 15OAf0. This model may be directly applicable to gravitational- wave detection of 
intermediate-mass black-hole mergers. 

PACS numbers: 04.25. Dm, 04.30. Db, 04.70.Bw, 04.80.Nn 95.30.Sf, 95.55.Ym 97.60.Lf 



I. INTRODUCTION 

Black-hole-binary mergers are a key target of ground- 
based and space-based gravitational-wave observations. 
The strongest radiation is produced just as the two black 
holes join to become one, and can only be fully under- 
stood through explicit numerical simulations. Since the 
first stable evolutions of black- hole-binary mergers [l|-lJl ; 
and after it was established that the gravitational wave- 
forms from these evolutions were universal, and consis- 
tent across codes and methodologies [5|-|7|, researchers 
have turned their attention to how the results of nu- 
merical relativity can most usefully be supplied to the 
gravitational-wave data-analysis community. 

After studying the equal-mass nonspinning case, re- 
searchers have had to address the complexity problem of 
more generic systems. Even allowing for simple scaling 
by total mass M = Mi + M2, and assuming zero eccen- 
tricity, such systems span a seven-dimensional parameter 
space: {77, Si, S2}, where 77 = MiM2/M^ is the symmet- 
ric mass ratio of the binary, and Si is the spin angular 
momentum vector of hole i. 

Early surveys of the waveform parameter space have 
restricted themselves to the 77-dependence of nonspinning 
systems. In Q, the authors investigated the multipolc 
structure of merger waveforms from such systems, not- 
ing that the strongest subdominant modes shared many 



characteristics with the dominant quadrupolc, and that 
they could be collectively described by an implicit ro- 
tating source model of the binary. The authors used this 
observation to construct a multi-mode gravitational-wave 
template family for such binary systems, as an alterna- 
tive to more usual effective-one-body (EOB) templates 

While we may assume that 77 and \Si\ remain essen- 
tially constant throughout inspiral and merger, the spin 
directions generally evolve, so a useful parametrization of 
the system should take care to distinguish components of 
the spin-direction space with physically distinct effects on 
the waveforms |11| . 

An obvious cut in parameter space to consider is that 
of spins aligned (or anti-aligned) with the orbital angular 
momentum. These systems will not precess, but exhibit 
observationally significant spin-orbit effects, distinguish- 
ing them from nonspinning binaries in their dynamics 
and resulting waveforms [l2| . High-accuracy waveforms 
from such evolutions have been produced and studied by 
several groups 13|- 16[ . Such systems have been partially 
characterized by [17j . using a variant of the frequency 
model from [8|. The frequency-domain phenomcnologi- 
cal templates of Ajith et al. have been extended to cover 
both mass ratio and total aligned spin [1^ [l^ , at least 
for the dominant modes, and attempts have been made 
to extend these to more generic systems [20, [2l| . 

A key result from our investigation of the dominant 



modes of nonspinning unequal-mass binary waveforms [g 
was that these modes had phases that evolved together 
in lock-step through inspiral and merger. This agreement 
was especially impressive for the £ = m modes, leading 
to the development of a heuristic picture of the binary 
system as a rigid rotator (at least in the adiabatic limit) 
driving gravitational-wave emission. We dubbed this the 
implicit rotating source (IRS) picture. 

A secondary result of this picture was the possibility 
of developing a simple model for the time-development of 
the dominant and leading subdominant modal frequen- 
cies UJim in terms of a single rotational frequency fl{t): 



^in 



mn{t). 



We also presented a simple model for the corresponding 
mode amplitudes, leading to the possibility of a new ap- 
proach to time-based gravitational-waveform templates. 
In fact, we developed such a template proposal, the 
IRS-EOB templates, as an alternative to the effective- 
one-body templates of [3, [lOl , which terminate the sig- 
nals by matching to a superposition of quasinormal-mode 
(QNM) frequencies. 

In this paper, we look at the dominant waveform modes 
from some aligned-spin systems, and ask the following 
general questions: Does the general IRS picture still 
hold? How do the features of aligned-spin mergers com- 
pare with those of nonspinning mergers? Can we quantify 
the main features of the merger ringdown with a simple 
analytic model? 

The rest of this paper is laid out as follows. In Sec. |lll 
we introduce the binary systems studied and the numer- 
ical methods used to simulate them. In Sec. IIIIl we 
present results for the final black-hole states and an IRS 
descriptive characterization of the waveforms from our 
numerical evolutions. In Sec. lIVl we analyze the late por- 
tions of these waveforms in more detail, and apply the 
analytic modeling approach of |8| to the dominant-mode 
frequencies and (with improvements) to the amplitude 
model, concluding with an explicit parametrization ap- 
proximating the (2,2) results of all our simulations. In 
Sec.|Vl we investigate the quality of the new models com- 
pared to the numerical waveforms in the context of the 
Advanced LIGO detector. We conclude with some dis- 
cussion in Sec. I VII Some extra detail on the convergence 
of the numerical simulations is given in the Appendix. 



II. SIMULATIONS 

To investigate the nature of aligned-spin binary wave- 
forms, we carried out a scries of numerical evolutions for 
equal-mass systems with zero spin (X1_00). spins aligned 
with the initial orbital angular momentum (X1_UU), anti- 
aligned (X1_DD), or mixed (X1_UD). We also re-ran, for 
purposes of comparison, the case of a 4:1 nonspinning 
binary (X4_00). 

The physical parameters of these evolutions are pre- 
sented in Table m The initial momenta of the equal-mass 



binaries, with the exception of Xl.UU, were chosen by 
integrating the post- Newtonian equations of motion, as 
outlined in [23, 01 , with spin contributions to the Hamil- 
tonian adapted from j24|-'28l| (although we work in the 
Arnowitt-Deser-Misner (ADM) gauge, the results from 
harmonic gauge using effective-field theory [26[ have been 
shown to be equivalent 123,1301), and the flux from l3ll |. 
For the X1_UU configuration, we used simpler quasicircu- 
lar initial parameters with no initial ingoing radial mo- 
mentum. For the X4_00 data, we retained the quasicir- 
cular initial parameters used in Q . 

The equal-mass runs were carried out with our Hah- 
NDOL evolution code [33 using the Paramesh mesh- 
refinement infrastructure [33^ . The new X4_00 data, how- 
ever, use Hahndol paired with the Einstein Toolkit [33| 
release of the Cactus Computational Toolkit [33 and the 
Carpet mesh- refinement driver |3a |. 

In all cases, the initial data are of the standard Brandt- 
Briigmann type |37| , using the Bowen-York l38| prescrip- 
tion for extrinsic curvature that exactly satisfies the mo- 
mentum constraint. We solve the remaining Hamiltonian 
constraint using the TwoPunctures spectral code [39J . 
This code also supplies the total ADM energy Madm of 
the system, as well as the individual "puncture ADM 
masses" Madm.i, to very high precision. We note, how- 
ever, that for highly spinning or boosted Bowen- York- 
type data, a measurable amount of radiation energy may 
be included in these puncture ADM masses, but then 
escape to infinity p, l4(]l - l42| : thus the initial puncture 
ADM mass may not be the optimal measure of pre- 
merger black-hole mass. These quantities are also listed 
in Table n 

To evolve these initial data, we employ the BSSNOK 
3-1-1 decomposition of Einstein's vacuum equations |43l - 
l45l . with the alternative conformal variable suggested in 
l46l448l | , constraint-damping terms s ugg ested in l49| , and 
the dissipation terms suggested in [50, [5i|. Our gauge 
conditions are the specific 1-t-log lapse and Gamma-driver 
shift described in [521 , which constitute a variant of the 
now-standard "moving punctures" approach [3|, [^ . 

The four equal-mass simulations - X1_00, X1_UU, 
X1_DD, and X1_UD - were conducted with the Hahn- 
dol/Paramesh version of our code using space-only 
adaptive mesh-refinement (AMR) with grids placed 
adaptively, based on curvature invariants [5[. The 4:1 
mass-ratio simulation X4_00 was carried out with the 
same evolution routines, now ported to run within Cac- 
tus/CARPET [3J-[30|, which applies mesh- refinement in 
time as well as in space. The initial grid structures for 
all runs are given in Table |IT1 For the equal- mass simula- 
tions, the highest-resolution regions closest to the punc- 
tures had a grid spacing of 3A^/160, M/64, or 3M/224 for 
the medium-, high-, and ultra-high- resolution evolutions 
(the ultra- high was performed for X1_UD only). For the 
X4_00 simulation, the grid resolution around the smaller 
hole was M/96, A'f/128, and M/160 for the medium-, 
high-, and ultra-high-resolution runs. 

The equal-mass simulations exhibit between second- 



TABLE I: Physical and numerical parameters of the initial data for all the runs presented, mi^p and m2,p are the bare puncture 
masses of the two pre-merger holes, rg and Po are the initial coordinate separation and (transverse) linear momentum, 
respectively, giving rise to a total initial orbital angular momentum Lq. Madm is the total energy of the initial data. The total 
infinite-separation total mass M of the system is estimated by Mah, the sum of the initial (apparent) horizon masses of the 
two holes, calculated at t = 100. finally, we quote the approximate observed eccentricity ([ij. 

run name 



mi^p = ?Ti,2,p 



Si, 



S2z 



ro Pof(xlO^) Po.(xlO^) 



Lq Madm Ej ^a: 



DM,i 



Mah en,, 



X1_00 


0.4872312 


0.0 


0.0 


11.0000 


9.00993 


X1_UU 


0.301805 


0.2 


0.2 


8.2013 


10.3248 


X1J3D 


0.390411 


-0.159125 


-0.159125 


11.9837 


8.83600 


X1_UD 


0.301805 


0.2 


-0.2 


11.0000 


9.00993 


X4_00 


0.7900, 0.1890 


0.0 


0.0 


8.4702 


6.95662 



7.09412 0.991092 0.990514 

0.0 0.846768 0.988459 

1.20000 1.058879 0.990453 

7.09412 0.991092 0.990024 

0.0 0.589240 0.992912 



1.000050 1.000050 0.002 

1.000908 1.000550 0.01 

0.998794 0.998686 0.01 

0.999222 0.998834 0.002 

1.000310 1.000315 0.02 



TABLE IL Initial grid structure of the different simulations. The leftmost number is the outer extent of the Cartesian grid, 
with resolution doubled (grid spacing halved) within each new refinement level. 



run name 



outer (fixed) grid structure inner (moving) grid structure 



X1_00, X1_DD, X1_UD 

X1_UU 

X4_00 (larger puncture) 

X4_00 (smaller puncture) 



[1536,768,384,192,144,72,24,12,8] [3.0,1.5,0.75] 

[1536,768,384,192,96,72,24,12,8] [3.0,1.5,0.75] 

[2048,1024,512,256,160,96] [20, 10, 5, 2.75, 1.5] 

[2048,1024,512,256,160,96] [20, 10, 5, 2.5, 1.25, 0.6875, 0.375] 



and fifth-order convergence for the Hamiltonian con- 
straint, while the momentum constraints only showed 
clear second-order convergence in the highest-resolution 
regions around the punctures. Nevertheless, waveform 
amplitudes and phases were sixth-order convergent over 
the majority of the evolution. The remaining simulation, 
X4_00, displays sixth-order convergence in waveform am- 
pHtude and phase until close to merger time. For details, 
we refer the reader to the Appendix. 

We use the AHFinderDirect code [H, HI] to locate 
the individual holes, as well as the final post-merger hole. 
We deduce the horizon mass from the horizon area Axr 
via Christodoulou's relation [55| 



A/Ih = A/f„. 



J2 



AMt 



where A/jn. = ^J Aj^y^/lQi: is the irreducible mass of the 
hole. We present the sum of the two horizon masses, 
A/ah = AfAH.i + A/ah,2 in Table HI and use it for time- 
scaling of gravitational waveforms. 

Following [23, llg , we estimate eccentricity using the 
variation in puncture orbital frequency ripu„c: 



en(t) 



V^ ^pu 



=) 



Zi Lr\r 



(1) 



where r^circ is a monotonic fit to ripunc, based on a simple 
post-Newtonian expansion. For a good fitting function, 
the residual en{t) should be a sinusoid of slowly decreas- 
ing amplitude, and period equal to the orbital period; 
the eccentricity is then the (nearly constant) amplitude, 
en, max- In practice, due in part to gauge-dependent be- 
havior in the puncture tracks, eo(i) is not perfectly si- 
nusoidal. Nevertheless, we quote the derived eccentricity 



measure for each run in Table ID This is higher than 
we would like for serious data-analysis applications, or 
for generating post-Newtonian-numerical-relativity hy- 
brid waveforms, and we could choose to reduce eccen- 
tricity through methods similar to those presented in [17[ . 
However, our primary purpose in this paper is to inves- 
tigate the bulk behavior of the waveform modes across 
configurations, and very low eccentricity does not appear 
to be necessary for this. 

To obtain gravitational waveforms from our simula- 
tions, we begin by calculating the "outgoing radiation" 
Weyl scalar ■04 |57|, corresponding to the tidal acceler- 



ations that are to be measured by gravitational wave 
instruments. ■04 is a complex quantity related to the 
wave strain h — h^ + ihy by two time-derivatives: 
■04 = — /i+ + i/ix [53|- We interpolate ■04 onto a set of 
coordinate spheres, and decompose the values on these 
spheres into spherical harmonics of spin- weight s = — 2, 



r'04(t,r,6',0) 



y ^Cjmjt, 



r)_2rr(0,0). 



em 



To obtain the harmonic modes of the strain h, there- 
fore, we must integrate Ci,m(t,r) twice in time, with in- 
tegration constants taken to yield zero strain long after 
the merger has taken place; we call this process "detrend- 
ing" the waveform. Currently, we use the Fourier-domain 
method of time- integration presented in [59| to produce a 
strain waveform h that is free of unwanted secular trends. 
This can also be written as a sum over modes: 



rft,(f, r, I 



E 

Im 



He,n{t,r)e'^'-^'^''^ -2Yr{0,<l>) . (2) 



In fact, we are most interested in an intermediate quan- 
tity, /i, which we call the strain-rate. This is of particular 
interest because it is most closely related to the rates of 
emission of g ravitational- wave energy and linear momen- 
tum [60|,|6l|. As with the strain and ip4, the strain-rate 
can be decomposed into spherical harmonics: 



r/i(t, r, ^, 0) = V A,,„(i, r)e'(v^'-^''^)+^/^) ^^y- 



E 



(3) 

where we have explicitly included a phase offset 7r/2 so 
that the remaining strain-rate phase ipim differs from the 
strain phase ^^,„ only by terms of 2.5PN order (see dis- 
cussion in jSTj). 

For the equal-mass cases, our extraction spheres were 
Text e {45M, 50M, 55M, 60M, 65M, 70M}, with a consis- 
tent extraction-region resolution of 6M/5, M, and 4M/5 
for central resolutions of 3M/160, M/64, and 3M/224, 
respectively. For the X4_00 case, the spheres were Text G 
{40M,50M,60Af,70A/,80M,90Af}. In these regions, 
the extraction-region resolution was M and 4M/5 for 
central resolutions of A//128 and A//160, respectively. 

In addition to errors in the strong-field region of the 
source, the extracted waveforms will also contain errors 
due both to finite extraction radius and finite grid resolu- 
tion in the extraction region. To mitigate the former, we 
have applied an extrapolation scheme to both waveform 
amplitude and waveform phase. Specifically, we assume 
a falloff model 



^rcxt — ^c 



a.2 



■, <Pr<,xt =(^00 + ^- (4) 

' cxt cxt 

Of all two-parameter Toxt-falloff models we have tried, 
this model gives the best fit to the amplitude and phase 
of the dominant (2, ±2) modes. Adding more terms to 
the falloff model will introduce overfitting errors, espe- 
cially given the limited range of our Text domain. This 
leading-order behavior is consistent with roxt-falloff pre- 
dictions of [62|. This model, however, seems inappro- 
priate for higher- frequency modes such as (4, ±4), where 
dissipation effects cause a general loss in amplitude. For 
these, we include the possibility of an additional term 
proportional to ^ext, at least for the amplitude: 



£\f 



a-ircxt + Aoo + 



0-2 
^2 



(5) 



We use this model for all modes with to > 3. We note, 
however, that diffusive effects should only act to decrease 
the amplitude. If a mode shows apparent growth that 
does not converge with some inverse power of rcxt, then 
it cannot be meaningfully extrapolated according to ([5]) . 



III. DESCRIPTIVE RADIATION 
CHARACTERIZATION 

As noted, the main objective of this paper is to char- 
acterize gravitational waveforms from aligned-spin merg- 
ers. In this section we present the main features of the 



radiation, in a spherical-harmonic modal decomposition. 
Our analysis follows the same approach developed in Q , 
which descriptively characterized the radiation from non- 
spinning mergers in terms of an implicit rotating source 
(IRS). In this approach, each modal waveform compo- 
nent is viewed as the trace of the dynamic development 
of one of a superposed set of source moments. To a very 
good approximation, each mode's radiation is circularly 
polarized, indicating rotational motion. This is registered 
in the waveform's rotational phase ^im = fim/'m'j while 
the modal amplitudes heuristically indicate the relative 
contributions of the source moments. 

Our goal is to build on the characterization of non- 
spinning mergers with additional details revealing the 
effects of aligned spins through the merger. We first 
characterize the raw content of the radiation in terms 
of energy and angular momentum, then comparatively 
examine how the modal amplitudes and rotation phases 
develop in time. 



Radiated Energy and Angular Momentum and 
Final States 



To calculate the rate of energy and angular momen- 
tum emission via gravitational radiation during merger 
and ringdown, we apply the following mode-summation 
formulas (see Appendix A of @): 



dE 

m 

dJ^ 
dt 



E 

Irn 



Ai 



167r 



Idtt 



"Pirn) : 



(6) 

(7) 



where we terminate the mode-sums at £ = 6 for the 
equal- mass cases, and £ ~ 5 for X4_00. We integrate 
the result in time to obtain the total energy -Brad and 
angular momentum Jz.md [x and y components arc zero 
by symmetry) emitted during the evolution. In princi- 
ple, these calculations could be performed with the full 
waveforms, rather than the {i, m) modes. In practice, 
however, we only output the mode-decomposed wave- 
forms for post-evolution analysis, and an accurate high- 
order time-integration of ipi within the evolution code is 
difficult. Additionally, post-evolution analysis with the 
waveform modes allows us to better control unphysical 
high-frequency noise. 

Rather than directly summing and integrating the Text- 
extrapolated strains and strain-rates, we instead inte- 
grate the finite- Text energy fluxes, and extrapolate the 
result according to the three-parameter fit 



AE, 



rad-Tc: 



AE, 



rad,c« 



£4 

„4 



Z^ + 4^. (8) 



and similarly for the radiated angular momentum, AJ^. 
However, when dissipation effects are significant, as for 
the X1_UU data, we must amend this assumption accord- 
ing to our model ([S]). Adding an Tcxt-proportional term 



TABLE III: Radiated energy and z angular momentum from 
all merger processes, in units of the infinite-separation total 
mass estimated by Mah (final column of Table |l]) . The pri- 
mary value is the r^xt-extrapolated value of the £ = 6 mode- 
sum of the integrals at highest spatial resolution, while the 
quoted uncertainty is the linear sum of three contributions: 
the standard error for the fit parameter for the £ = 6 mode- 
sum at the highest resolution; the difference between the £ — 4 
and £ — 6 mode-sums at this resolution; the difference be- 
tween the £ — 6 mode-sum result for the highest and next- 
highest resolutions (for X4_00, the £ = 5 mode-sum was used 
instead of ^ = 6). 



run name A_Erad(A/AH) 



AJz,Tiid{Mj^ii) 



X1_00 0.038547 ± 0.000244 0.367786 ± 0.001117 

X1_UU 0.075636 ± 0.001413 0.482665 ± 0.003874 

X1_DD 0.027240 ± 0.000219 0.293735 ± 0.001084 

X1_UD 0.039792 ± 0.000440 0.373655 ± 0.001581 

X4_00 0.014437 ± 0.000104 0.136347 ± 0.000853 



to the strain-rate amplitude will introduce several new 
terms to a quadratic- in-amplitude quantity like E'rad- 
However, since this many terms are impossible to fit cred- 
ibly with only six extraction radii, we instead extrapolate 
the waveform modes first according to ^ (for m < 4) 
and ([5]) (for rn. > 4), and then perform a mode-sum of 
the resuh. The results are given in Table Hill 



Now we present our estimates of the final state of the 
post-merger Kerr holes, encoded in the two parameters 
A'/f and a = Sz/Mf. Our estimates are derived from a 
number of sources, and are tabulated in Table ITVl 

Most directly, the columns marked Aff^rad and arad are 
derived from simple conservation of energy and angular 
momentum: 



Mf,,ad 
arad 



Jo — AJ2 rad 

-"^faad 
Lo + Siz + S'22 — A J^.r 



ad 



^"f,rad 



(9) 



(10) 



where A£',ad and AJ^ j-ad are taken from Table Hill and 
the remaining quantities are as in Table HI 

We can compare with an end-state model based on fits 
to a range of numerical mergers. One such model for final 
mass, appropriate for equal-mass systems, was given by 

A^f,AEi/MAH = l-po-Pi(ai + a2)-P2(ai+Q;2)^, (11) 



where a a = I'S'a/M^J is the initial dimensionless spin of 
hole A, and the fitting parameters are (again, determined 
by comparison with numerical data): 



Po = 0.04826 ±0.00027, pi 
P2 = 0.00485 ±0.00025. 



0.01559 ±0.00026, 



(12) 



We note that the uncertainties on the parameters are 
incomplete, with an undetermined (but presumably neg- 
ligible) post-Newtonian component. 

For the final spin, one model with just enough com- 
plexity for our data sets here was given by [65l |66| ^: 

"AEI = a + Sirja'^ + s^ri'^a + tQria + 2V3ri + t2'ri'^ +^7]^, 



q-^ai + a2 



q- 



1 



(13) 



where the coefficients {54,55,^01^21^3} were determined 
by comparison with numerical data: 

Si = -0.1229 ± 0.0075, S5 = 0.4537 ±0.1463, 
io = -2.8904 ± 0.0359, t2 = -3.5171 ±0.1210, 
h = 2.5763 ±0.4833. (14) 



In Table IIVI we present final masses and spins derived 
from values derived from Eqs. (|lip and ([T^ . with un- 
certainties due only to the parameter uncertainties in 
Eqs. (dSj) and dH]). Unfortunately since Eq. ^ only 
applies to equal-mass cases, we cannot use it to estimate 
the X4_00 end-state mass. The more complicated formula 
found in Lousto et al. |64| covers more generic binaries. 



but with larger uncertainties. 

We note also from Table HVl that the two cases X1_DD 
and X4_00 have the same final spin, within the quoted 
uncertainties. Thus we might expect similarities in the 
ringdown portion of their waveforms across all impor- 
tant modes, though the extent to which each quasinor- 
mal mode (QNM) is excited will be different in the two 
cases. 



B. Multipolar amplitudes 

In Ref. H ■^fi found strong similarity in the peak-scaled 
modal amplitude development through the peak for a 
range of nonspinning mergers over a range of masses, and 
somewhat rougher similarity among the different modes. 
For nonspinning mergers, the dominant modes were gen- 
erally those with £ = m, and these modes were neatly 
described with the IRS heuristic. 

Strain-rate amplitudes for the strongest modes of our 
new simulations are shown in Fig. [Tj For all equal-mass 



^ Other models for the post-merger mass arc available; see, for 
instance, Tichy & Marronetti [6311 and Lousto et al. |64| . 



^ Note that we have adapted Eq. 4 of [65l | to match our convention 
for q. 



TABLE IV: End-state Kerr parameters (M, a) of post-merger holes. Mf^rad and ftrad, and associated uncertainties, are derived 
from radiation balance (|9I10|I - see Table Hill Aff.AH and OfAH come from the AHFinderDirect code [53. l54| and the Hahndol 
spin calculator [671 ]: quoted uncertainties are a combination of the post-merger variability of the irreducible mass and spin and 
the difference between the measured mass and spin from the highest and second-highest resolutions. Aff,AEi(AfAH) and oaei 
use the numerically tuned formulas (|13p and pip due to [14l. [6a. [6a| : quoted uncertainties here are due to uncertainty in the 
fitting coefficients HH), (fT4)l . 



run name 


Mt,rad(MAH) 


Qrad 


A/f,AH(MAH) 


Qah 


A/f,AEl(A/AH) 


aAEi 


X1_00 


0.9519 ± 0.0002 


0.6878 ± 0.0013 


0.95165 ± 0.00001 


0.68644 ± 0.00001 


0.9517 ± 0.0003 


0.68646 ± 0.00004 


X1_UU 


0.9123 ± 0.0014 


0.9165 ± 0.0055 


0.91164 ± 0.00013 


0.90720 ± 0.00015 


0.9144 ± 0.0008 


0.9114 ± 0.0264 


X1_DD 


0.9645 ± 0.0002 


0.4825 ± 0.0012 


0.96303 ± 0.00002 


0.48140 ± 0.00012 


0.9637 ± 0.0006 


0.4794 ± 0.0256 


X1_UD 


0.9514 ± 0.0004 


0.6847 ± 0.0019 


0.94996 ± 0.00001 


0.68408 ± 0.00002 


0.9517 ± 0.0003 


0.6865 ± 0.0243 


X4_00 


0.9782 ± 0.0001 


0.4726 ± 0.0009 


N/A 


N/A 


N/A 


0.4748 ± 0.0093 



simulations, the strongest subdominant modes are (4,4) 
and (3,2); other modes never attain 0.1% of the (2,2) 
power (equivalently, 3% of the (2,2) amplitude). For 
X1_DD, the (2,2) mode is even more strongly dominant: in 
this case, all subdominant modes other than (3,2) show 
significant power only at very late times. At i? = 45Af 
the (2,0) mode shows an amplitude similar to the weaker 
of the modes shown here, but this is sensitive to the ex- 
traction radius (our procedure for detrending the strain- 
rate doesn't work well for m = 0). 

It is worthwhile to briefly consider how the modal com- 
position varies with aligned spin and mass ratio, as shown 
here. Note that the (4,4) mode amplitude is roughly the 
same for all cases shown here, varying even less than the 
(2,2) amplitudes, as was already seen for the nonspinning 
runs investigated in Q. For the equal- mass cases, the 
(3,2) mode amplitude roughly equals the (4,4) mode at 
peak, but for aligned (anti-aligned) spins it is enhanced 
(suppressed) approaching the peak. For the 4:1 mass- 
ratio X4_00, the odd £ = m make significant contribu- 
tions unseen for equal masses, more so than in the X1_UD 
asymmetric-spin case. For asymmetric mergers of either 
kind, £ = |r?T,| -f 1 modes are also significant. 

Figure [5] overlays the (2, 2) (top panel) and (4, 4) (bot- 
tom panel) amplitude peaks of all cases (suitably time- 
shifted and rescaled) to compare their relative sharpness. 
It is interesting that the "down-down" peaks of X1_DD are 
narrower than the "up-up peaks" of X1_UU. The steeper 
slope on the i > side can be tied to the generally faster 
fall-off in QNM modes for the prograde modes of the 
much more slowly rotating black hole generated by the 
down-down merger. The peaks remain roughly symmet- 
ric, with a faster rise as well. This is particularly striking 
for the subdominant modes. 



C. Waveform phasing 

In our studies of nonspinning mergers [8|, we found 
strong correspondence in phase development among the 
different modes, interpreted as near- "corotation" of the 
implicit-source moments. Specifically, all significant 
modes displayed a common rotational phase $£,„ up to 



the time of peak power at merger, deviating by less than 
0.025 rad during that time. After the merger, the de- 
viations between modes increased, but for the i = m 
modes this deviation was very slow, ^ 1 rad over the first 
lOOA'f following merger. For the weaker i ^ m modes, 
the phasing began to differ somewhat earlier and was in 
some cases less cleanly described by the IRS heuristic. 

Modal phase comparisons are more challenging for 
ec^ual-mass spinning mergers than they were for nonspin- 
ning unequal-mass mergers. As noted above, the symme- 
tries of the configuration and the weakness of radiative 
spin effects in the inspiral yield only a few significant 
modes and even these tend to be weak, subject to com- 
petition with noise in the simulations, and likely more 
sensitive to subtleties in the choice of spherical-harmonic 
basis. 

The left panel of Fig. |3] shows the phasing of several 
strain-rate modes for the up-up case X1_UU. Generally, as 
was the case for nonspinning systems Q, the different 
(£, m) modes remain approximately in phase up to the 
merger. As with the nonspinning mergers, the i = m 
cases show the best agreement for t < 0. In this case 
though, because exchange symmetry excludes the odd- 
m modes, we only have two i = m modes to compare 
up to £ = 5. Two modes present - (4,2) and (5,4) - 
have amplitudes below our 3% cut-off in Fig. [1] Such 
small amplitudes introduce a lot of noise in the mode's 
phase; we include the phase only when it begins to show 
acceptable continuity. 

We take a closer look at the relative phasing in the 
right panel of Fig. [3l where we present the difference be- 
tween each of the three strongest subdominant modes - 
(4,4), (3,2), and (5,4) - with the dominant (2,2) mode. 
Generally the phase differences decrease going from a fi- 
nite extraction radius (dashed curves) to Text — ^ oo (solid 
curves). 

Looking at the inspiral portion (t < 0) of the phase 
first, we see that the (2,2) and (4,4) rotational phases 
agree within ~ 0.05 rad, with a marginal improvement 
when we extrapolate r^xt -^ oo. For the (3,2) mode, 
there is a roughly constant offset of about 0.15 rad af- 
ter Text-extrapolation. Due to the short extent and noisy 
nature of the reliable (5,4) mode phase, it is difficult 



g 0.01 




^ 



0.1 



0.01 



0.001 



0.0001 



jit. l ■ ^f^L^^^nf 




-500 



peak ^ AH' 



FIG. 1: Multipolar strain-rate amplitudes for X1_UD (top 
panel), X1_UU (second), X1_DD (third), and X4_00 (bottom 
panel) , evaluated at finite extraction radii (45M for the equal- 
mass, 50M for the X4_00 case). The lack of symmetry for 
X1_UD and X4_00 means additional excited (odd-m,) modes. 
For the X1_UD case, we overlay circles to show the nearly iden- 
tical even-TTi mode amplitudes of X1_00. In each panel, we 
omit subdominant amplitudes that never exceed 3% of the 
(2,2) mode. 




FIG. 2: Comparison of "sharpness" of amplitude peaks for 
(2, 2) (top panel) and (4, 4) (bottom panel) modes of all cases. 
Each amplitude has been rescaled by its maximum value. 



to extract a definite phase offset; it appears to be in the 
range ^ 0.05 — 0.10 rad after extrapolation. However, the 
phase offset has also flipped sign during extrapolation in 
Text, indicating that we may not know the correct phase 
to high accuracy. This may not be resolved simply by 
increasing grid resolution; we have seen similar extrap- 
olation sign-flips for the (5,4) phase offset in our X1_UD 
simulations, even at the "ultra-high" 3M/224 resolution. 
As with the nonspinning case, the i ^ m modes show the 
largest offset from the (2, 2) modes, and are most affected 
by Text -extrapolation effects. 

We note that post-Newtonian theory predicts for 
nearly constant phase offsets between modes during late 
inspiral; these come in at 1.5PN order for certain modes 
(see, for example, the polarization amplitudes given in 
[69|). However, they are small compared to the phase 
differences shown here - less than ~ 0.03 rad up to lOOM 
before peak. 

Looking now at the post-merger period {t > 0), the 
phase agreement remains quite tight, better than that 
seen in the nonspinning mergers. In the IRS interpre- 
tation, all modes in this case remain nearly rotationally 
locked right through the merger. Note in particular that 
the phase difference between the (2, 2) and (3, 2) modes 
is roughly constant for t > 0. The (4, 4) and (5, 4) modes 
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FIG. 3: Left: Rotational phase $(f) from puncture tracks and different multipolar strain-rate components extrapolated to 
Text — > oo for the up-up case X1_UU. The weakest two modes, (5,4) and (4,2), are only measurable for about 200M before 
merger. Right: The difference in phase with the (2,2) mode for the next three most important modes only: (4,4) (circles), 
(3, 2) (squares), and (5, 4) (diamonds). In each case, we show the difference at Text = 45M (dashed lines), and r^xt — ;■ oo (solid 
lines) . 



are also in phase with each other at Text = 45Af ; they 
develop a phase offset when extrapolated to r^xt ^^ o^i 
but maintain the same slope. All modes agree within ^ 1 
rad even 60M after peak. 

To understand this tight phase agreement, we may look 
to perturbation theory for the post-merger Kerr hole. For 
rapidly spinning black holes, QNM frequencies depend 
primarily on m, approaching Mwqnm = Tfijl in the a — )■ 
A'/f limit [ZOilZll I which suggests a tighter coupling for the 
modes in this case. However, the final spin of the post- 
merger Kerr hole for X1_UU, a « 0.91, is not close enough 
to this extremal limit to explain the phase agreement we 
see. 

In the left panel of Fig. |4] we show the phasing of several 
modes for the down-down simulation X1_DD. Again, the 
(4, 2) and (5, 4) modes are weak, and yield reliable phases 
only from ^ lOOM before peak. Looking at the right 
panel, the (4,4) mode extrapolated phase difference is 
^ 0.02 rad, while there is a small offset of about 0.1 rad 
between the (3, 2) and (2, 2) modes. As with the X1_UU 
mode, extrapolation in Text appears to increase the phase 
offset. There is a slight drift among the i = m modes, 
similar to that seen in other cases. Indeed, the phase 
difference between (2.2) and (4,4) is nearly identical to 
that of the X1_UU case. The £ ^ m modes this time 
show significantly varying frequencies (slopes). This is 
unsurprising, consistent with the differences among the 
leading normal QNM mode frequencies. 



IV. MODELING THE NEAR-MERGER 
WAVEFORM 

In this section we undertake a more quantitative study 
of the late-time waveforms through merger and ringdown, 
following the general approach in [Sj for an explicit quan- 



titative representation of the frequency development. We 
extend the previous work on modeling the amplitude with 
additional parameters to allow more precise fits at rela- 
tively early times. 

The waveform phasing examined in the last section is 
fairly featureless. The phase is monotonie, slowly devel- 
oping curves with a gentle elbow at merger. This sim- 
plicity is a result of the slow secular development of the 
underlying circular motion which generates the radiation. 
It also suggests that we may quantify the phase develop- 
ment with just a few parameters. 

Following the approach in [Sj we probe more deeply 
into the phasing by taking a time-derivative to study the 
frequency evolution. Common features are found among 
the leading waveform modes and across a range of merg- 
ers, allowing the results to be summarized with a simple 
parametrization. With the same general frequency model 
as in Q for nonspinning mergers we can also describe the 
phasing of spinning black-hole mergers. 

In Fig. [21 we compare the dominant-mode frequencies 
of the three equal-mass cases presented above. Since 
odd-m modes are suppressed by symmetry, the rel- 
evant modes are the (2,2) (quadrupole) , (4,4), and 
(3, 2) modes. Unsurprisingly, the frequencies are consis- 
tently higher throughout merger for more-aligned spins, 
with the final plateau value matching the dominant 
quasinormal-mode (QNM) frequency. For the (3, 2) 
modes, there is significant deviation from the smooth 
frequency development generally expected according to 
our IRS heuristic; this amounts to a large bump in the 
frequency during the plateau phase. Similar deviations 
in the (3,2) modes have been noted previously [^, \V\ . 
Such effects may arise through mode mixing with the 
(2, 2) mode [73], which could arise through ambiguity in 
the shape of the sphere on which the radiation is mea- 
sured, or on the use of (spin- weighted) spherical harmon- 



50 



40 



fe 



30- 



20- 





' 


'1 


- 


■ puncture track 

- (2,2) 

- (4,4) 


A 




- (3,2) 

- (5,4) 


(/// 


' 


- (4,2) 


]//// 


_,^'- 



-200 



-100 



0.5 






-0.5 



<^^(4,4)-(2,2)(45M) 
--(3,2)-(2,2)(45Af) 
^^(5,4)-(2,2)(45M) 
.„ (4,4)-(2,2) (extrap.) 
(3,2)-(2,2) (extrap.) 
<^^ (5,4)-(2,2) (extrap.) 



-500 



'"U'^^ah) 




FIG. 4: Left: Rotational phase "l>(i) from puncture tracks and different multipolar strain-rate components extrapolated to 
Text — !■ oo for the down-down case X1_DD. The weakest two modes, (5,4) and (4,2), arise very suddenly at late times, and 
are only measurable for about lOOM before merger. Right: The difference in phase with the (2,2) mode for the next three 
most important modes only: (4,4) (circles), (3,2) (squares), and (5,4) (diamonds). In each case, we show the difference at 
''cxt = 45A'f (dashed lines), and rcxt — ;■ oo (solid lines). 



ics, rather than the spheroidal harmonics appropriate for 
the perturbation theory in which the QNM frequencies 
are defined 
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The precise cause and mechanisms of 
this mixing are open questions, which we hope to return 
to in future work. 

Also included in Fig. [5] are the equivalent frequen- 
cies for the 4:1 nonspinning merger X4_00 (note, how- 
ever, that that merger had significant odd-m modes not 
present in the equal- mass cases here) . The spins of the 
anti-aligned X1_DD initial data were chosen to yield the 
same final Kerr parameters (mass, spin) as the X4_00 
data, according to ((T^ . As the Kerr parameters deter- 
mine the QNM frequencies of each mode, it is not sur- 
prising that the X1_DD and X4_00 frequencies level off to 
the same value after merger. What is interesting is the 
difference in behavior approaching this final state. For 
t < -20A/, X4_00 hews closely to the nonspinning X1_00. 
The latter could be expected given the similarity in phas- 
ing upon approach to merger for nonspinning mergers |8[ . 
At the latest times the frequency development is deter- 
mined primarily by the parameters of the final black hole 
formed, while additional parameters become important 
as we look back to earlier times. 



A. Modeling the Rotational Frequency 



In [8| , we introduced the following empirical model for 
the rotational frequency ^im = 'J-'fm/'Ti in a short time- 



window around the merger: 

i7(i) = r!f(i-/(i)) 



f{K,b,to;t) = ^fl + - 

Z \ K 



{1+k) 



(l + ie-2(t-to)/6 
V 1^ 



(15) 



,(16) 



where the dimensionless parameter c = Clob/^i replaces 
the "maximum frequency slope" Oq corresponding to the 
peak chirp rate. Within this general framework, there are 
up to five free parameters for the frequency development: 
c, K, 6, to J and fif. 

In our previous investigations [8|, this functional form 
worked well in fitting the dominant frequencies of a se- 
quence of nonspinning binaries with mass ratios in the 
range {1.0,6.0}. Unsurprisingly, fif was found to be 
consistent with the quasinormal frequency of the post- 
merger Kerr hole. More interestingly, b was also found to 
be approximately consistent with the quasinormal damp- 
ing time, meaning that at late times the frequency ap- 
proaches its limiting value exponentially at the same rate 
as amplitude squared. It was also found that the dimen- 
sionless ratio AffSlo/^f ~ 0.021 across all cases. 

We usually apply PS|) as an "orbital frequency" , which 
is scaled from the gravitational- wave frequency by the az- 
imuthal mode number m. The formula may be applied, 
with similar results, to strain, strain-rate or -04 wave- 
forms. 

We consider three increasingly constrained classes of 
fits of this form. The most general is a free fit for 
all five parameters. Second, we test the late-time fre- 
quency/amplitude relationship noted in [8| with a fit 
where b is constrained to agree with the late-time ampli- 
tude fall-off rate (and thus with the QNM fall-off rate). 
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FIG. 5: Waveform frequency Loira for the dominant modes of 
the equal-mass evolutions. The upper panel shows the dom- 
inant (2, 2) mode, while the middle and bottom panels show 
the next strongest modes - (4,4) and (3,2), respectively. We 
also show the corresponding frequencies for the 4:1 nonspin- 
ning merger X4_00. At early times, this tracks the X1_00 wave- 
form, while during merger it approaches the final frequency 
of the X1_DD case. 



Finally we suggest a general fit by which all parameters 
(other than to) are derived from the final black- hole mass 
and spin. 

We focus on the strain-rate fit, as it provides a good 
compromise between the base-line drift error that affects 
the strain waveforms, and the higher level of noise in the 
■04 waveforms. Figure [S] shows the result of this pro- 
cedure for the Text-extrapolated waveforms of the X1_UU 
(top) and X1_DD (bottom) runs (note that the numerical 
data have been down-sampled by factors of 20 or more for 
clarity of presentation). At the level of precision acces- 
sible by eye, all fits appear nearly perfect after t > —20. 
Parameter fits conducted only over times t > — 20 typi- 
cally do not extrapolate well to earlier times; a fit over a 
wider range, extending over t > — 40 appears to be very 
good over this entire region though there is some slight 
degradation in the quality of the fit near < t < 10. For 
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FIG. 6: Rotational frequency fits for the (2, 2) modes of the 
X1_UU (top) and X1_DD (bottom) runs, extrapolated to Text -^ 
oo. The different fits differ primarily in the early part of the 
comparison window around ipcak, which runs from — 20Mah 
to -|-40Mah for standard fits, and from — 40Mah to -|-60AfAH 
for wide fits. 



the dashed curves labeled "free fit", all five parameters 
in (|T6l) were fit freely against the numerical data. In the 
curves labeled "& fixed" , we test the hypothesis that the 
exponential decay of frequency evolution is related to the 
amplitude fall-off rate, fixing b first by a fit to the expo- 
nential decay rate in the mode amplitude data, before 
fitting the other parameters according to the frequency 
data. Though the frequency data are then fit against four 
parameters instead of five, the result is still a very good 
fit, justifying the assumption. The amplitude/frequency 
relationship is discussed more in the next section. With 
b based on the amplitude data, we record the best-fit val- 
ues of the parameters c, k, b, J7f, and <o in Table [V] The 
final curves in Fig. |S] (labeled "constrained") only fit tg, 
with all other parameters pre-set, as discussed below in 
Sec. UVCl 



B. Amplitude Modeling 

Following [8|, our strategy is to describe the wave am- 
plitudes in relation to the frequency. This is loosely mo- 
tivated by the idea that frequency evolves in response to 
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loss of energy and angular momentum, but energy and 
angular momentum fluxes are dependent on the wave am- 
plitude. For nonspinning systems, we previously found 
that dJ/dfl was slowly varying in the merger, and could 
be approximated as a constant. The result was enough to 
provide a coarse quantitative description of the merger- 
ringdown amplitude development in reference to the fre- 
quency development. Here we extend that model, in- 
troducing additional parameters to allow more precise 
quantitative description of the numerical results. 

Now assume the waveform strain-rate amplitude takes 
the general form; 



Airnif)^AoimPif)\ f (t) , 



(17) 



where the adjusting function P{x) is some function that 
approaches unity as x — >■ 0. Then assuming that the 
strain amplitude Hgm satisfies |iifo„| ^ \Hemip'lj^\, which 
is true at all points of the numerical waveform, we can 
find an approximate expression linking strain and strain- 
rate: 



H-hn^tm- 



(18) 



Combining this with (|17p . we can model the strain am- 
plitude as: 



(mflf)-i(l - ./)-' Ao,™ P{f)J fit) 



Hofm{l-f)-'P{f)\ fit) 



(19) 



Similarly, the amplitude of the tJj^ i£, m) mode would be 
modeled as: 



Cimif)^Coimil-f)Pif)\ fit) 



(20) 



where we are still using the strain sign convention for 
phasing of modes (that is, positive-m modes have positive 
frequencies). 

Parameters including Aq and other parameters in the 
definition of Pif) allow some tuning for various cases 
considered here while preserving the general approach in 
0. Concretely, consider 



PifY 



N 

E 

n=l 



a, 



(/'" 



?2n+2 



(21) 



The simplest possibility, with A^ = yielding P = 1, 
would imply that dE/dw is constant, i.e. that the sys- 
tem loses radiative energy in linear proportion to the late- 
time frequency decay to the quasinormal-ringdown rate; 
this is close to the amplitude model used in Q. Equa- 
tion (j2ip is consistent with quasinormal ringing radiation 
in the / — > limit (assuming frequency model parameter 



b = I/Qwqnm) and can be adjusted for deviations earlier 
in the waveform where 1 > / > 0. The restriction to 
even powers was motivated by an empirical observation 
that the first helpful correction seems to be at second 
order, and the resummed powers in the summand, yields 
a more generally regular result as / ^- 1. Going back to 
times more than 20M before merger the model becomes 
unrealistic. The model amplitude begins small at early 
times, growing exponentially toward the peak. 

In practice, we find that we get a good approxima- 
tion for the merger-ringdown part of the radiation in the 
(2,2) modes by keeping one term in the expansion (|2ip . 
and fitting for Aq and ai. Then using / from ()16p . the 
complete amplitude model used is 






\rh,,if)\' ^ Al 



fit) 



l + «i /2-/ 



(22) 



The result of this procedure is shown in Fig. [71 for the 
Toxt-extrapolated waveforms of the X1_UU and X1_DD runs. 
For the less-constrained fits, b was first determined using 
data in a window from 20A'/ah to SOMah! this value was 
then fixed, and data over the larger window — 20A/ah to 
IIOMah were used to determine Aq and ai. Though not 
shown in the figure, we see somewhat less accurate fits 
with the model for the amplitudes of the X4_00 case, with 
differences before and near peak at the ^ 5% level. 



C. Constraining the Models 

Working with the results of the free fits for frequency 
and amplitude, we note the approximate constancy of 
the parameters c and k across all cases. Additionally, we 
note that the final frequency and decay parameters Of 
and b are close to the expected QNM values. Thus we 
may be able to reduce considerably the number of free 
parameters needed for the models. 

Similar to the dimensionless scaling used for the c 
parameter above we seek a scaling of ai in terms of 
the QNM "quality factor" Q = 3fiwQNM/29a;QNM, as a 
simple dimensionless number dependent on the spin of 
the final hole. We find that the results for our equal- 
mass cases roughly scale with Q^, with the mean result 
ai ~ 72.3/g2. 

We now perform a more constrained version of the fits, 
fixing the parameters b, flf to their QNM values, and 
replacing c and k with their average values from Table FVl 
and setting the frequency parameters and ai as outlined 
above. The combined set of constrained parameters is: 



c = 0.252 , K = 0.426 , r^f == 5Ra;QNM/2 
b = 1.0/9cjQNM,ai = 72.3/Q2. 



(23) 



Thus we arc left with just two free parameters to fit: to 
and Aq. These (as well as the constrained parameters) 
are recorded in Table IVII 
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FIG. 7: Amplitude fits for tfie (2, 2) modes (extrapolated to 
roxt — s> oo) of the X1_UU (top) and X1_DD (bottom) runs. The 
two fits in each panel differ in whether they fit the parameter 
ai or just Aq. The fit windows here emphasize the late tail 
of the amplitude, running from — 20AfAH to +120AjrAH- 
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FIG. 8: Frequency (top panel) and amplitude (bottom panel) 
fits for the (4, 4) mode (extrapolated to Text — >■ oo) of the 
X1_DD run. Again, the numerical data are indicated by circles, 
with a free fit over the parameters represented by the dashed 
line. 



We may consider attempting to constrain the remain- 



ing parameters as well. The Aq parameter has units Mj 



1/2 



and seems to scale approximately with ilj 
mean fit for the equal-mass cases we get 



1/2 



9.9r]il 



1/2 



Using the 



(24) 



We include rj in the fit since the overall amplitude coeffi- 
cient must vanish linearly as the mass-ratio goes to zero. 
Though we have not focused on mass-ratio dependence, 
this scaling is consistent with the nonspinning 4:1 result. 



ting analysis, the (4, 4) mode is plicnomcnologically sim- 
ilar to the (2, 2) mode. Figure |8] shows the frequency 
and amplitude fits for the (4,4) mode of the X1_DD run. 
The frequency fit is clearly still very close to the numeri- 
cal data over the domain of interest, but the amplitude's 
overall peak is ^ 10% too low, with a poor fit to the 
slope of the numerical data before the peak. This sug- 
gests that our ansatz for the mode amplitude does not 
carry over to subdominant modes, and requires further 
work. Nevertheless, we will see in the next section that 
the dominant mode may already be useful in detection 
studies. 



D. Subdominant modes 

As noted above the most significant modes for equal- 
mass mergers are the (4, 4) and (3, 2) modes. Even these 
have amplitudes of no more than about one-tenth that of 
the (2, 2) mode. While the (3, 2) mode shows more com- 
plicated features that do not lend themselves to this fit- 



To contrast the quality of the fit performance for the 
(2, 2) and (4, 4) modes, we present in Fig.[S]the associated 
strain-rate (real parts) over the range of the fit. We also 
plot in the upper panel the (2, 2) mode resulting from a 
fully constrained model for amplitude and phase, using 
Eqs. dMl-dllD. 
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TABLE V: Best-fit values for the unconstrained parameters c, k, b, Q[, and to for tiie frequency model (|15|) . and of Ao and 
ai for the amplitude model (|17[1 . All fits are over a time window from tpcak — 40M to tpcak + 60M. Quoted uncertainties are 
the direct sum of three terms: uncertainties in the highest-resolution fits; differences between best-fit values for rext — >■ oo and 
''oxt ~ 45 Af (40M for X4_00); differences between best-fit values at highest and next-highest resolutions runs. 



run name c 



to 



Qi 



Ao 



ai 



X1_00 0.2489 ± 0.0040 0.421 ± 0.015 -3.77 ± 0.17 11.685 ± 0.025 0.27655 ± 0.00021 1.270 ± 0.022 6.64 ± 0.35 

X1_UU 0.2500 ± 0.0021 0.401 ± 0.009 -2.15 ± 0.59 14.296 ± 0.030 0.37317 ± 0.00070 1.627 ± 0.064 2.46 ± 0.70 

X1_DD 0.2626 ± 0.0062 0.473 ± 0.027 -4.90 ± 0.12 11.203 ± 0.034 0.23805 ± 0.00093 1.157 ± 0.009 11.12 ± 0.12 

X1_UD 0.2458 ± 0.0038 0.407 ± 0.014 -3.57 ± 0.15 11.634 ± 0.024 0.27643 ± 0.00017 1.263 ± 0.016 6.27 ± 0.38 

X4_00 0.2343 ± 0.0033 0.439 ± 0.009 -4.73 ± 0.52 11.381 ± 0.033 0.23380 ± 0.00075 0.740 ± 0.017 9.96 ± 1.26 



TABLE VI: Values for the parameters c, k, b, fif, and to for the frequency model (|15|) . and of Ao and ai for the amplitude 
model (|17|l . Unlike in Table [V] only to and Ao are freely fit; the remaining parameters have been fixed, as given in Eq. (|23|) . 
All fits are over a time window from tpoak — 40M to tpeak + 60M. Quoted uncertainties are the direct sum of three terms: 
uncertainties in the highest-resolution fits; differences between best-fit values for r-oxt — > oo and roxt = 45M (40M for X4_00); 
differences between best-fit values at highest and next-highest resolutions runs. 



run name c 



to 



nt 



Ao 



Ct-1 



X1_00 0.252 0.426 -3.99 ± 0.29 11.712 0.27661 1.271 ± 0.013 6.7934 

X1_UU 0.252 0.426 -2.31 ± 0.58 14.404 0.37133 1.633 ± 0.032 2.4924 

X1_DD 0.252 0.426 -4.41 ± 0.36 11.222 0.23820 1.143 ± 0.006 9.9784 

X1_UD 0.252 0.426 -3.91 ± 0.30 11.681 0.27661 1.272 ± 0.006 6.8296 

X4_00 0.252 0.426 -6.67 ± 0.72 11.465 0.23305 0.732 ± 0.009 9.9867 



V. FAITHFULNESS OF THE FREQUENCY 
MODEL 



One way to quantify how much of the merger informa- 
tion we have captured by the modehng above is to com- 
pare the results of the model with the original fully nu- 
merical waveforms in a detector context. We consider an 
explicit waveform model restricted to the (2, ±2) modes, 
which contain most of the power. The waveform phase is 
derived from integrating the model IRS frequency given 
by (jTS]) and (fTB]). while the amplitude is given by ([22]) . 

In total there are seven parameters in these expres- 
sions, and one additional parameter ipo arises as an in- 
tegration constant in deriving the phase from our fre- 
quency model. Drawing on the results of Sec. IIVCI five 
of these parameters {c, k, b, flf, ai} are specified by (^5)) 
as functions of the final black hole's leading quasinormal- 
mode frequency, thus reducing these free parameters to 
functions of the final black hole's mass and spin. Beyond 
these the only remaining parameters are {ipo^to-, Aq}^ cor- 
responding to phase and time references, and an overall 
amplitude scale. For the XIJDD data, we can sec the re- 
sulting waveform in the top panel of Fig. [9] For compar- 
ison we also show the results of a "free" fit, where the 
frequency was fit to the NR data without the constraints 
in (|23p . and the amplitude was fit based on those results 
without constraint on ai. This is not a completely free 
fit, since the frequency parameters are fit without regard 
for the consequences on the amplitude. 

We calculate the mismatch for the Advanced LIGO 
detector [73. "Mismatch" here is defined as the devia- 



tion of the normalized overlap integral from unity, usually 
optimized over free parameters such as overall phase and 
arrival time iTal: 



mismatch = 1 — max ■ 



{h„i{Xi)\he 



A. ^{h„,{X^)\K^i\^)){he\K 



(25) 



where the frequency-space inner product (-I-) is the over- 
lap between two signals, defined as [Tg 



{hi\h2 



hifMir + hi{fyh2U) 



= 4Re 



^n(/) 



■df 



Snif) 



df 



(26) 



where hi(f) and /i2(/) arc the Fourier transforms of the 
signals, and Sn{f) is the (one-sided) noise power spectral 
density of the detector we're interested in; in this case, 
we take the ideal form Sn{f) for Advanced LIGO given 
in Appendix A of [Tj]: 



Snif) = So{x 

+ 111 ( 1 



-4-14 _ 52^-2 



X 

~2 



X 

~2 



.(27) 



where x = f / f^, Sq = lO'^^ and /o = 215 Hz. 

For this test, we pick black- hole binaries of total mass 
M e {40Afo,300MQ}, at a fixed distance of IGpc from 
the detector. Furthermore, we assume the system is ob- 
served along the polar axis. The resulting mismatch with 
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FIG. 9: Real part of the (2, 2) (top panel) and (4, 4) (bottom 
panel) strain-rate modes for the X1_DD run (extrapolated to 
Text — >■ CX3). Again, the numerical data are indicated by circles, 
with a free fit over the parameters represented by the dashed 
line. For the (2, 2) mode, we also include the result of the 
fully restricted waveform (see Eqs. ((23)l ) as a continuous line. 



the quadrupole NR signal is shovifn for all the runs (X1_00, 
X1_UU, X1_DD, and X4_00) in Fig.[10l Common to ah cases 
is a sharp fallofF (that is, improvement) in mismatch as 
the system mass increases: from > 25% at M = AOMq 
to < 1% for M > 20QMq. This trend is expected: over- 
all physical frequencies scale inversely with system mass, 
so while the last few pre-merger orbits' worth of radia- 
tion for a 40 Mq system might fall in Advanced LIGO's 
most sensitive frequency band, only the higher-frequency 
merger and ringdown might lie in the same band for a 
200Mq system. 

In Fig. [TT] we again show this mismatch for the X1_00 
configuration, along with three other mismatches: with 
the full NR signal (all modes) along the polar axis, and on 
the equatorial plane, and also against the unconstrained 
model, using "free" -fit parameters (Table ly]). This last 
mismatch has been scaled up by a factor of 100 for visi- 
bility. 

For these results we have compared with just the (2,2) 
component of our numerical simulation waveforms. If ap- 
plied in an actual observation there would be additional 
power, perhaps at the level of up to several percent in 
other harmonics. This model makes no attempt to fit 



rt 0.2 

s 



1 


o-exi 


UD 


q-eXI 


UU 


*-*Xl 


DD 


^"^X4 


00 




M{Mq) 

FIG. 10: Mismatch (f25|l between quadrupole (2, ±2) con- 
strained IRS and numerical-relativity waveforms for all simu- 
lations in the context of the Advanced LIGO detector, where 
the system was observed along the polar ax;is at a distance of 
IGpc. 
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FIG. 11: Mismatch ^5^ for the X1_00 configuration, be- 
tween quadrupole (2, ±2) constrained IRS waveforms and (a) 
quadrupole numerical-relativity (NR) waveforms along the 
polar axis (circles), (b) full NR waveforms along the polar 
axis (squares), (c) full NR waveforms on the equator (stars), 
and (d) the unconstrained IRS waveforms along the polar axis 
(triangles). All mismatches were calculated in the context of 
the Advanced LIGO detector, at a distance of IGpc. 



those contributions. With a little work we could extend 
our constrained model to approximate the contributions 
of these other modes, but this would necessarily require 
dependence on many additional parameters, including in- 
formation about the component masses and spins and the 
relative orientation of source and detector. 

We note that the the "sweet spot" of the Advanced 
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LIGO sensitivity curve (|27|) is around 200 Hz. For sys- 
tems at the low- mass end of our plot, 4OAf0, this fre- 
quency range is accurately handled by post-Newtonian- 
based waveforms. Thus a full waveform model appropri- 
ate for such low-mass systems should really be a combi- 
nation of our merger-ringdown model with a PN-based 
inspiral. By neglecting this here, and integrating over 
the full band, we will suffer from junk numerical radia- 
tion and windowing artifacts at the lower-mass end. It 
is reasonable to suspect that we would get a consider- 
ably lower mismatch if we restricted our integration to 
frequencies above Mu)22 ~ 0.15. The estimates presented 
in Figs. fTOUTT] are therefore likely to be conservative. 

With one detector, signals based on the waveform 
model wc have constructed here depend only on the in- 
trinsic parameters describing the final black hole, phase 
and time references and overall amplitude, just five of the 
17 parameters describing generic black- hole merger ob- 
servations. Our results suggest that dominant features 
of the powerful merger-ringdown radiation may be de- 
scribed with little or no reference to the details of the 
component black holes. This provides a complementary 
description of the merger to those based on the binary 
inspiral parameters. Even without a long inspiral lead- 
in, such models may be useful in detecting gravitational 
waves from high-mass mergers. Our waveforms are based 
on the same parameters as those in ringdown-based ap- 
proaches to merger observations in ground-based detector 
pipelines J78l48l| and may be useful in future versions of 
these searches. 



VI. DISCUSSION 

In this paper, we have investigated the mode- 
decomposed gravitational waveforms resulting from the 
merger of aligned-spin black-hole binaries. Our primary 
purpose was to establish how well the implicit rotating 
source (IRS) picture of the binary as a GW source - first 
suggested in Q in the context of nonspinning holes of 
comparable mass ratios - holds in this different branch 
of parameter space. 

Based on these investigations, we note that the modal 
structure of aligned-spin mergers is like that of the 
equal-mass nonspinning configuration, dominated by the 
(2, ±2), (3, ±2), and (4, ±4) modes. These modes still 
display IRS-type behavior, featuring common rotational 
phase evolution with little offset through late-inspiral, 
merger, and into ringdown. The peak modal amplitudes 
are similar to those for the nonspinning case, though the 
duration of the peak region (which was roughly indepen- 
dent of mass ratio) is extended for aligned spins (and 
shortened for anti-aligned spins). A similar timescale de- 
pendence is seen in the rise to peak frequency. 

In applying our late-merger frequency model (with a 
slightly modified paramctrization) to these new cases, we 
have found that the model still performs well for the dom- 
inant modes. We enhance the original mode amplitude 



model of [8| to achieve improved behavior, at least for 
the leading (2, ±2) modes; however it yields up to '^ 10% 
mismatches at and before peak for the next most impor- 
tant modes. 

For the (2, ±2) modes at least, we have attempted to 
constrain all parameters explicitly with reference only 
to the state of the final black hole (i.e., its dominant 
quasinormal-mode complex frequency). With these con- 
straints, wc have reduced the additional free parameters 
to just three: ioj the time of maximum chirp-rate, <poj 
the phase offset, and ^o, the mode's amplitude scale. 
This description provides an approximate fit to the late 
part of the waveforms for all our simulations, including 
equal- mass spinning cases, and the 4:1 nonspinning case. 
Moving back in time to earlier points before the merger, 
the quality of this fit degrades and other physical details 
of the premerger binary become more significant. We 
see evidence of this when comparing the X4_00 (4:1 non- 
spinning) and the X1_DD (equal-mass down-down spins) 
configurations, which result in the same final spin. 

We have quantified the quality of this approximation 
by calculating Advanced LIGO fitting factors. For sys- 
tem masses of ^ 15OA/0, wc have found mismatches of 
< 5% between the full numerical-relativity waveform and 
the (2, ±2)-mode-only model waveforms. 

Our results suggest that an approach to gravitational- 
wave observation templates with parameters tied first 
to the structure of the final black hole may be useful 
for Advanced LIGO observations of intermediate-mass 
mergers. These would be an alternative to the time- 
domain effective-one-body (FOB) templates of |82h84{ . 
and the frequency-domain "phcnomenological" templates 
of [1^, [l^ , similar to ringdown searches currently being 
applied to LIGO data J78l - l8l| . In future work we plan 
to investigate the quality of this model, or its extensions, 
for a broader variety of mergers, including processing and 
eccentric configurations. 

In this explicit modeling, we have focused on the dom- 
inant (2, 2) mode. A similar model incorporating full 
multi-mode information can also be applied to com- 
plete inspiral-merger-ringdown waveform templates, as 
was done for nonspinning systems in [8|. Future work 
on this topic will focus on improving and extending the 
amplitude model to cover multiple significant modes of 
the merging binary. 



Acknowledgments 



This work was supported by NASA grants 08-ATFP08- 
0126 and 09-ATP09-0136. Resources supporting this 
work were provided by the NASA High-End Comput- 
ing (HEC) Program through the NASA Advanced Su- 
pcrcomputing (NAS) Division at Ames Research Center. 



16 



Appendix: Convergence 

In this appendix, we present the convergence proper- 
ties of the evolution fields and extracted waveforms. Our 
presentation will model that of Q . 

The simulations X1_00, X1_DD, and X1_UD used identical 
numerical methods and grid structures, with finest (near- 
puncture) resolutions of 3M/128, 3M/160, and Af/64, 
and wave-extraction- region resolutions of 3M/2, 6M/5, 
and M, respectively. The X1_UD simulation was also car- 
ried out at an ultra-high resolution of 3Af/224 (wave- 
extraction resolution 6A//7). 

The remaining equal- mass simulation, X1_UU, uses 
identical numerical methods at the same resolutions, but 
had a different grid structure in the wave-extraction zone. 
This has been seen to result in higher noise levels in wave- 
form quantities, but should not affect the overall conver- 
gence properties. We will use the X1_UD resolution to 
assess convergence levels for the Hahndol/Paramesh 
code. 



1. Constraints 

To establish constraint convergence in the equal-mass 
runs, we look at the Ll-norms of the Hamiltonian and 
momentum constraints both in the strong-field region 
and in the region where the waveforms are extracted. 

Figure [T^ shows the Ll-norm of the Hamiltonian con- 
straint. In the upper panel, level 13 (the region just out- 
side the punctures) demonstrates between fourth- and 
fifth-order convergence; in the lower panel, level 8 (con- 
taining the wave-extraction spheres) demonstrates di- 
minished convergence - between second and third order. 
Note that the resolution in these outer regions is much 
lower than in the crucial high-resolution regions, where 
the black holes reside. It is our understanding that er- 
rors in these distant regions are dominated by the effects 
of uninteresting short-wavelength features (particularly 
gauge modes), which propagate out from the center and 
become poorly resolved in the coarse regions. 

Convergence is more difficult to establish for the mo- 
mentum constraint. Figure [T51 shows the behavior of the 
X component's Ll-norm in level 14 (the region contain- 
ing the punctures), and level 8 (containing the wave- 
extraction spheres). The very inner zone (level 14) - 
displays clean behavior: all components of the momen- 
tum constraint are roughly 2.6-order convergent. Once 
we move outside this finest region, however, convergence 
drops down precipitously. We note, however, that the 
momentum constraint in general is two to three orders of 
magnitude lower than the Hamiltonian constraint, which 
suggests that there is a small amount of low-order er- 
ror present in all constraints, but which is dominated 
by higher-amplitude (but convergent) error only in the 
Hamiltonian constraint. 

■Gonstraint violation information was not available for 
the Cactus-based X4_00 simulations. 




t(M) 



FIG. 12: Convergence of the Hamiltonian constraint's Ll- 
norm for the X1_UD simulation: between fourth- and fifth- 
order convergence in level 13 (upper panel); between second- 
and third-order convergence in level 8 (lower panel). 







' 1 ' 




0.00015 




med 

f(2)*high 

f(3)«high 


- 




1 RvAA/viAAAAAaaa I 




WmmiiM 


0.0001 - 


vvvvvwwi/w! 

1 1 1 1 1 1 1 1 1 1 1 






' 1 








med 




" 




f(0)*high 

f(l)*high 


" 


2 0.5 - 


\MHi 




T| l^l\ jiMiMlmii 




vmB^h 




n m^ jsi>mKiil~ 


L 


1 1 1 1 1 1 1 1 1 1 1 1 1 



FIG. 13: Convergence of the momentum constraint's Ll-norm 
for the X1_UD simulation: between second- and third-order 
convergence in level 14 (upper panel) ; without any established 
convergence order in level 8 (lower panel). 
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FIG. 14: Amplitude convergence for the (2, 2) (upper panel) 
and (4, 4) (lower panel) modes of the X1_UD simulations, based 
on the central resolutions 3Af/160, M/64, and 3Af/224. In 
each case, the (3Af/160 — 3Af/224) difference has been scaled 
up assuming fourth- and sixth-order convergence. Though 
noisy, the amplitude differences appear to be consistent 
with sixth-order convergence throughout the inspiral, until 
^ 150Af before merger. 



FIG. 16: Amplitude convergence for the (2, 2) (upper panel) 
and (4, 4) (lower panel) modes of the X4_00 simulations, based 
on the central resolutions Af/96, Af/128, and A//160. In each 
case, the (Af/96 — A4"/160) difference has been scaled up as- 
suming fourth- and sixth-order convergence. The amplitude 
differences appear consistent with sixth-order convergence un- 
til approximately 6QM before peak, when the rate declines to 
close to fourth-order. 
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FIG. 15: Phase convergence for the (2, 2) (upper panel) and 
(4, 4) (lower panel) modes of the X1_UD simulations, based 
on the central resolutions 3A//160, Af/64, and 3A//224. In 
each case, the (A//64 — 3Af/224) difference has been scaled up 
assuming fourth- and sixth-order convergence. The phase dif- 
ferences appear to be consistent with sixth-order convergence 
throughout the inspiral, merger, and ringdown. 
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FIG. 17: Phase convergence for the (2, 2) (upper panel) and 
(4, 4) (lower panel) modes of the X1_UD simulations, based 
on the central resolutions Af/96, A///128, and A'//160. In 
each case, the (A4"/96 — A'f/160) difference has been scaled 
up assuming fourth- and sixth-order convergence. The phase 
differences appear to be consistent with sixth-order conver- 
gence through late inspiral, declining to closer to fifth-order 
at merger and ringdown. 



In Figs. [13] and [121 we demonstrate sixth-order con- 
vergence for the (2,2) and (4,4) modes' amplitudes and 
phases, respectively. 

For the 4:1 data (Fig. fTH]) . we see generally cleaner 
waveforms, but also a large oscillation in errors until 
about 400A^ into the evolution. After this, conver- 
gence appears to be sixth-order until close to ampli- 
tude peak time, when it declines to fourth-order. This 
may be because the overall error is dominated by un- 
certainties in the merger time, which is determined by 



the fourth-order-accurate Runge-Kutta time-integration 
scheme. Unfortunately, we could not disentangle this ef- 
fect to sufficient accuracy to establish sixth-order conver- 
gence through the peak time. 

The 4:1 waveform phase evolution (Fig. fTT)) also seems 
to display sixth-order convergence until the merger, when 
it declines to fifth-order. We note that the scale of the 
errors is generally less than half those of the X1_UD data 
from Fig. [T5] 
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